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Abstract 

In the present paper, we derive a closed-form solution of the multi-period portfolio 
choice problem for a quadratic utility function with and without a riskless asset. All 
results are derived under weak conditions on the asset returns. No assumption on the 
correlation structure between different time points is needed and no assumption on the 
distribution is imposed. All expressions are presented in terms of the conditional mean 
vectors and the conditional covariance matrices. 

If the multivariate process of the asset returns is independent it is shown that in the 
case without a riskless asset the solution is presented as a sequence of optimal portfo- 
lio weights obtained by solving the single-period Markowitz optimization problem. The 
process dynamics are included only in the shape parameter of the utility function. If a 
riskless asset is present then the multi-period optimal portfolio weights are proportional to 
the single-period solutions multiplied by time-varying constants which are depending on 
the process dynamics. Remarkably, in the case of a portfolio selection with the tangency 
portfolio the multi-period solution coincides with the sequence of the simple-period solu- 
tions. Finally, we compare the suggested strategies with existing multi-period portfolio 
allocation methods for real data. 
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1 Introduction 



Nowadays, the portfolio selection problem plays an important role in financial research. A 
number of papers are devoted to questions like, e.g., how an optimal portfolio can constructed, 
monitored, and/or estimated by using historical data (see, e.g., Alexander and Baptista (2004) , 
Golosnoy and Schmid (2007), Bodnar (2009)), what is the influence of parameter uncertainty on 
the portfolio performance (cf., Okhrin and Schmid (2006) , Bodnar and Schmid (2008)), how 
do the asset returns influence the portfoHo choice (see, e.g., Jondeau and Rockinger (2006), 
Mencia and Sentana (2009), Adcock (2009), Harvey et al. (2010), Amenguala and Sentana 
(2010)), how is it possible to estimate the characteristics of the distribution of the asset returns 
(see, e.g., Jorion (1986), Wang (2005), Frahm and Memmel (2010)), how can the structure of 
optimal portfolio be statistically justified (Gibbons et al. (1989), Britten- Jones (1999), Bodnar 
and Schmid (2009)). 

In a seminal paper from 1952 H. Markowitz presented the idea of an optimal portfolio 
selection by taking into account the trade-off between the portfolio expected return and its risk 
which is measured by the portfolio variance. The idea of Markowitz 's approach is to minimize 
the portfolio variance for a given level of the expected return. This method is equivalent to the 
so-called mean- variance utility maximization problem. Although the suggested methodology 
is quite simple, it provides us the most commonly used solution of the single-period (static) 
portfoho choice problem that remains very popular today (see, e.g. Brandt (2010) ). 

Although the case of a long-term investment horizon is of greater importance in practice 
much less has been done in that area. The first formulation of the multi-period portfolio 
selection problem has already been given in the book of Markowitz (1959) followed by the 
papers of Mossin (1968), Samuelson (1969), Merton and Samuelson (1974). Although it is 
heavily discussed in recent literature (see, e.g., Li and Ng (2000), Steinbach (2001), Leippold, 
Trojani and Vanini (2004), Brandt and Santa Clara (2006), Qelikyurt and Ozekici (2007), Skaf 
and Boyd (2009)), to the best of our knowledge, a closed- form solution is not available in the 
general case up to now. Only under the assumption of independence closed-form solutions 
are presented by Li and Ng (2000) and Leippold et al. (2004). For more general models, the 
solution is frequently determined by a numerical procedure (see, e.g. van Binsbergen and Brandt 
(2007)). Brandt and Santa Clara (2006) suggested a solution of the multi-period portfoho 
selection problem assuming that the portfolio weights can be presented as a linear function of 
certain state variables. This assumption leads to a simplification of the optimization problem. 
Note that the solution is only a local maximum which could differ from the global one. Finally, 
the solution of the multi-period portfolio selection problem in continuous time is given by Duffie 
and Richardson (1991), Basak and Chabakauri (2010), Ai't-Sahalia et al. (2009) among others. 
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We contribute to the existing literature by deriving the closed-form solution of the dynamic 
portfolio choice problem with and without a riskless asset under rather weak assumptions. 
The only conditions imposed on the distributions of the asset returns are the existence of the 
conditional mean vectors and of the conditional covariance matrices. No assumptions about the 
correlation structure between different time points or about the distribution of the asset returns, 
like normality, are needed. The suggested method can be applied for both stationary and non- 
stationary stochastic models. The results are obtained assuming that the investor makes his 
decision on the basis of the quadratic utility function. This is one of the most commonly used 
procedures since the paper of Tobin (1958) where it is shown that the BernouUi principle is 
satisfied for the mean- variance solution only if one of the following two conditions is valid: the 
asset returns are normally distributed, which is rarely the case in application, or the utility 
function is quadratic. On the other hand, the quadratic utility function is usually considered 
as a good approximation of the other utility functions (cf. Brandt et al. (2006)). Moreover, 
under the additional assumption of independence we show that at each time point the optimal 
multi-period portfolio weights can be presented in a similar way as the optimal single-period 
portfolio weights. Both representations differ in the coefficient of risk aversion. Finally, if the 
asset allocation is based on the tangency portfolio we prove that the solution of the multi-period 
portfolio selection problem is the same as that obtained by solving the single-period problem 
at each time point. 

In an empirical study we apply the obtained results to real data by comparing the perfor- 
mance of the suggested strategies with existing multi-period portfolio allocation methods. It is 
shown that the multi-period portfolio strategies based on the approximative solutions perform 
very well for different values of the coefficient of risk aversion and different investment periods. 

The rest of the paper is organized as follows. In Section 2, we introduce the multi-period 
portfolio choice problem for a quadratic utility function without a riskless asset. The main 
results of this section are given in Theorem 2.1 (Section 2) where a closed- form solution of 
the optimal portfolio weights is given. In Corollary 2.2 we present the results of Theorem 
2.1 in terms of the classical one-period Markowitz's solution for independent observations. In 
Section 3, the solution of the multi-period portfolio selection problem for a quadratic utility 
function with a riskless asset is presented (see Theorem 3.1). As in the case without a riskless 
asset the solution for independent returns is proportional to the solution of the corresponding 
single-period problem at each time point. The process dynamics exclusively influence the 
proportionality constant (Corollary 3.2). A very interesting result is obtained for the tangency 
portfolio. In Theorem 3.3 we prove for independent returns that the multi-period portfolio 
selection problem for the tangency portfolio is equivalent to a sequence of the single-period 
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problems. The same solution is obtained in both cases. Because the tangency portfoho is, 
usually, considered as a market portfoho (see, e.g. Britten- Jones (1999)) in capital asset pricing 
theory the obtained result is of great importance for practitioners. The findings of the empirical 
studies are presented in Section 4. The paper concludes in Section 5. All proofs are given in 
the appendix (Section 6). 

2 Multi-Period Portfolio Choice Problem for a Quadratic 
Utility Function: Without Riskless Asset 

In this section we derive a closed-form solution of the multi-period portfolio choice problem with 
k risky assets for the case of a quadratic utility function. Let = (-^t,i, -^t,2, • • • , -^t,fc)' denote 
the vector of the returns of k risky assets and let E(Kt\J^t-i) = and Var(Kt\J^t-i) = S^. Tt 
denotes the information set available at time t. is assumed to be a positive definite matrix. 
Note that for deriving the closed-form solution of the multi-period portfolio selection problem 
in the present case we need neither any assumption on the correlation structure between dif- 
ferent time points nor any distributional assumption. It is only demanded that the conditional 
covariance matrix of the asset returns exists. The solution exclusively depends on the condi- 
tional mean vector and the conditional covariance matrix. These quantities can be calculated 
depending on the underlying model of the asset return process. For instance, if the investor 
assumes that the asset returns follow a vector ARMA-GARCH process then the expressions for 
the conditional mean vector and the conditional covariance matrix can be directly obtained by 
applying the well-developed theory of the multivariate autoregressive processes and the mul- 
tivariate conditionally heteroscedastic autoregressive processes (see, e.g., Brockwell and Davis 
(1991), Engle (1982, 2002), Bollerslev et al. (1988), Bauwens et al. (2006)). 

Let Wt = {wt^i^ Wt^2^ • • • , '^t,/c)' denote the vector of portfolio weights in period t. Note that 
w^l = 1 where the vector 1 denotes a /c-dimensional vector whose components are all equal to 
1. Then the wealth of the investor at time t is given by 

Wt = Wt.,{i + w;_iX,) = Wt.,wi,Xt , (1) 

where = 1 + X^. Then E(Kt\J^t-i) ^ P^t ^ ^ + l^t Var(Kt\Tt-i) = ^t- Later on, 
we make use of X^ instead of X^ for simplifying the presentation of the obtained results. In 
this section, we deal with an investor who invests his money exclusively into k risky assets and 
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whose investment strategy is based on the quadratic utihty function given by 



U{W,) = W,- -Wf , (2) 

where a > is the slope parameter of the quadratic utihty function. The corresponding relative 
risk aversion coefficient (RRA) is given by 

It = -V7T (3) 

1 — aWt 

specifies the attitude of the investor toward risk. Brandt and Santa Clara (2006) considered jt 
as a constant obtained by substituting = 1 in ([3]). We follow this procedure in the empirical 
part of the paper by choosing the value of a in the definition of the utility function ^ in such 
a way that 7 G {5, 10, 15,20}. 

The portfolio is allocated at time point and thereafter reconstructed at time 1, . . . , T — 1. 
The planning horizon contains T periods. The investor gets his reward after the final period at 
time T. The aim of the investor is to maximize his final utility of wealth at time point T, i.a. 

y(0,H^o,^o) = max Eo[U{Wt)] . (4) 

Here Et{') stands for the expectation given the information set J-'t available at time t. We 
assume that the short-selling is allowed, i.e. could be negative. This problem can be solved 
recursively. 
Let 

V{t,Wt,J't)= max EtpiWr)] (5) 

and let denote the optimal portfolio weights at time point T — t + 1 which is in general 

a function of optimal weights of the next periods, i.e. of w^_^^2 5 ^T-t+s^ ^t-i- Following 
Brandt and Santa-Clara (2006) and Pennacchi (2008) the optimization problem ^ can be 
solved by applying the following Bellman equation at time point T — t 



V{T-t,WT-t,J'T-t) 



max Ex-t 

WT-t:w^_(l=l 

max Ex-t 

WT-t:Wy_(l=l 



max 



Ex-t+iWiWx)] 



V{T-t + 1, W^r-tW^-i+iXT-t+i, 



(6) 
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for t = 1, . . . , r. Note that 



V{T,Wt,J^t) = U{Wt). 



(7) 



First, we derive the closed-form solution for the weights at period T — 1 and then formulate 
the general solution in Theorem 2.1. At time point T — 1 the application of the terminal 
condition ([t]) leads to 



V{T, Wt, Tt) = U{Wt) = Wt- ^W^ 



(8) 



and, hence, 



V[T - 1, Wr-i, J^T-i) = max Et-y 

{wT-i:w^_-^l = l} 



max Et-\ 

{wt-i:w^_-l1=1} 



a 



max 

{wt-i:w^_-l1=1} 

max 

{wT-i:w^_-^l=l} 



/ ^ ( ^ ~ \ 2 1 



where At = St + fiTi^T 

The first order conditions for the period T — 1 are 



(9) 



Wt-i^It oW^-^At^t-i + A1 = and l^w^-i 1 



(10) 



where A is a Lagrange multiplier. Solving (10) with respect to A and w^-i leads to the portfolio 
weights for period T — 1 given by 



w 



'IV A- ' 



T-l 



— 7— i 1 — TTr — QtAt with Qt = Arj} ^ ^ 1 



In Theorem 2.1, the expressions of the optimal portfolio weights w^_^ at periods to T — 1 
are presented. 

Theorem 2.1. Let = (X^^i, X^^2, • • • , -^r,/c)^ r = 0, . . . , 6e a random return vector of k 
risky assets and let £'(X^|J>_i) = /x^ and ]/ar(X^ |J>_i) = where Z5 positive definite. 
Then for all t = 1, . . . ,T t/2,e optimal multi-period portfolio weights for the period T — t are 
given by 



(12) 



6 



with 



Q 



\T-t+l 



1 

■T-t+l 



1M"^ 1 



(13) 



St + AtAt /^^ t = 1 

^T-t[VT-i+2XT-mX^_,+i] for t = 2,...,T, 



(14) 



.At for t = l 

^T-t[i?r-t+2XT-i+i] /or t = 2,...,T 

and 

RT-,v.= ^-^%if^. l/.-« = — ^ /or t = 2,^^^,r. (16) 

The proof of the theorem is given in the appendix. Theorem 2.1 shows that the structure 
of the optimal portfoho weights remains the same during the investment period. 
It is important to note that if the terminal condition ([7| is replaced by 

V{T, Wt, J^t) = U{Wt) = Wt- -{Wt - ET-i{WT)y 

then its solution is the same as in Theorem 2.1 except the expression at time point T — 1. Here, 
we get 

11 ^1 1 1 

An important property of the optimal weights is shown in Corollary 2.2 that is proved in 
the appendix as well. Here, it is shown that under the assumption of independence the optimal 
weights have the same structure as the ones that are obtained by solving the single-period 
portfolio selection problem at each time point with the common quadratic utility function 
given by 

Q/ 

max /xjw — — w'StW subject to w'l 1 . (17) 



7 



In the following we use the notations 

l^SrV, 1 



GMV,i 



Corollary 2.2. Let = (X^^i,X^^2, • • • ,-^r,/c)^ r = 0, . . . 6e a sequence of the indepen- 
dently distributed vectors of k risky assets and let E(X.r) = jn^ VarCX^r) = Sr where S^- 
is positive definite. Then for all t = 1, . . . , T the optimal multi-period portfolio weights for the 
period T — t are given by 

y-l -I 

w^-t = T + a^_^+iQr-f+iMT-f+i ^ (18) 



and 



^T-t+i - ^T-t+i 1,^-1 — ^ — U^; 



T 

n ) ~ 1 ~ RcMV.T-t+l 



= ^^^^ (20) 



ai = 



1 + Rgmva 



(1 + RgMV^iY + (1 + '5i)VbMy,i 



The results of Corollary 2.2 are similar to those of Li and Ng (2000) who solved the multi- 
period portfolio-selection problem assuming that the asset returns are independent. Corollary 
2.2 states that for solving the multi-period portfolio selection problem with the quadratic utility 
function ^ it is enough to solve the single-period problem given in (17) at each time point 



tG{0,l,...,T — 1} and then to make an adjustment in the expression for the coefficient of the 
investor's risk aversion a. It is very remarkable that the optimal weights at time T — t only 
depend on the mean vector and the covariance matrix at time T — t + 1 and the risk aversion 
aT-t^i' The mean vectors and the covariance matrices at time points r > T — t + 1 have an 
influence on the optimal weights only over the quantity aT-t+i- 



Because in the case of the single-period portfolio selection problem the solution of (17) lies 
on the eflicient frontier where a specifles the location of the optimal portfolio within the eflicient 
frontier, the same behavior can be observed in the case of the multi-period portfolio allocation 
problem. The only difference is that the efficient frontier is time- varying since the mean vectors 
[1^ and the covariance matrices are used for its construction. Then the optimal portfolio 
is obtained by choosing as specifled in Corollary 2.2. In the special case when the process 
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{X^-} consists of independent observations with time-invariant mean vector and time-invariant 
covariance matrix the efficient frontier remains the same during the whole investment period. 
Nevertheless, the optimal portfolio obtained by solving the multi-period portfolio selection 
problem is time- varying because is a function of future wealths and it is not constant. 

Hence, in the case of investing exclusively into risky assets and the solution of the multi- 
period portfolio selection problem for the quadratic utility function is not equivalent to the 
solution of the corresponding T single-period allocation problems. Even if we assume that 
{X^} consists of independent observations and the mean vector and the covariance matrix are 
constant during the investment period, the risk aversion is time-varying. 

The results of Corollary 2.2 possess another important application. Using the relationship 



between the mean- variance utility optimization problem (17) and the Makowitz optimization 
problem we get the formulation of the later in the multi-period case. At time point T — t it is 
given by 

min v^'YiT-t+i^ subject to (21) 

n ^i] — RcMv.T-t+i ) ^T-t+l 



1 



1 + ST-t^l 



w'l = 1 
for t = 1,...,T. 

3 Multi-Period Portfolio Choice Problem for a Quadratic 
Utility Function: With Riskless Asset 

In this section we present a closed-form solution of the multi-period portfolio selection problem 
with k risky assets and one riskless asset. Let be the random return vector containing of k 
risky assets as defined in Section 2 and let r/^t be the return of the riskless asset at time t. Let 
Wt = {wt^i^ Wt^2^ • • • , '^t,/c)' denote the vector of portfolio weights in period t invested into the k 
risky assets. Then the evolution of the investor's wealth is expressed as 

Wt = Wt-1 (1 + Tf^t + w;_i(Xi - rf,tl)) = Wt-1 [Ru + wj.iiti) , (22) 

where i?/^t = 1 + and = X^ — The conditional mean vector of X^ is /i^ = /x^ — 

and its conditional covariance matrix is given by S^. In this section we consider an investor 
who invests into k risky assets and one riskless asset with the investment strategy based on the 
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quadratic utility function given in The optimization problem is given by 



y(0,Wo,J-o)= max Et[U{WT)] 

Let 



(23) 



V{t,Wt,J't)= max EtpiWr)]. 



(24) 



The problem (23) is solved recursively by applying the following Bellman equation at time 
point T — t 



V{T-t,WT-t,J'T-t) 



= max Et-i 
= max Et-i 

WT-t 



max ET-t^i[U{WT) 
V{T-t + 1, Wr-t (vf^T-t + w;.^,+iXT-t+i ) , J^T-t+i) 



(25) 



for t = 1, . . . , T. Note that for t = the equality ([T]) holds. 

As in Section 2, we first derive the closed-form solution for the weights at period T — 1. After 
that the general solution is derived (see Theorem 3.1). At time point T — 1 the application of 
([T]) leads to 



y (T, Wt, Tt) = U{Wt) = Wt- ^W^ . 



(26) 



Let At = + (j^tI^t Because 



Et-1 (WT-l{Rf,T + w^_;l^t)^ = VaiT-i (^VFt-i(^/,t + w^_;l^t) 



we get 



F(T-1,Wt-i,J^t-i) = max Et-i Et[Wt - -W} 

{wt-i} L 2 



max £^T-i 



max 

Wt-1 



= max 

Wx"-! 



Wt-1 {Rf,T + ^'t-iI^t) - ^^T-l (w^_;l^tWt-1 + i^/,T + 2i?/,TW^_i/i^) 
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The first order condition for the period T — 1 is 



Wr-ifJ^T - oiW^_^ (^At^Wt-i + Rf^Ti^r) = . (27) 
This leads to the following expression of the portfolio weights for the period T — 1 

(rfib 

In Theorem 3.1, the optimal multi-period portfoho weights are presented. The proof of the 
theorem is given in the appendix. 

Theorem 3.1. Let X^- = (X^-^i, • • • , -^r,/c)'; ^ — 0, ...,T^ he a random return vector 
of k risky assets^ let rf^r be the return of the riskless asset, and let E{lLr\Tr-i) — fjir 
Var(Kr\J-'r-i) — where is positive definite. Then for all t — 1,...,T the optimal 
multi-period portfolio weights for period T — t are given by 

^T-t = \^^rrJ, ^ ( n ^f^' ) ~ Rf,T-t+i ^T-t+iMr-t+i (29) 



with 



, St + UtI^t for t — 1 ^ ^ 

AT-t^r={ ^ ^^^^ . . ^ , (30) 

ET-t[{l - 5T-t+2)XT-t+iX^_,^J for t = 2, . . . , T 



t^T t i={ '^^ ^"""^ (31) 

^ Er-t ((1 - ST-t+2)^T-t+i) for t = 2, . . . , T 



and 



ST-t+2 — /ir-t+2^r-t+2Mr-t+2 /^^ t — 2, . . . , T . (32) 

As in the case without a riskless asset the expression of the optimal portfolio weights for 
each period T — t looks like the solution of the single-period portfolio selection problem at time 
point T — t. In Corollary 3.2, a stronger result under the assumption of independence is proved, 
namely that the obtained weights are proportional to the weights obtained by solving 

max /xjw — — w'S^w . (33) 



11 



Corollary 3.2. Let = (X^^i,X^^2, • • • ,-^r,/c)^ r = 0, . . . be a sequence of the indepen- 
dently distributed vectors of k risky assets^ let rf^^ be the return of the riskless assets and let 
E(X.r) = l-ir l/ar(X^) = where is positive definite. Then for all t = 1^ . . . the 
optimal multi-period portfolio weights for period T — t are given by 



w 



T-t ~ 



(34) 



where 



T 

aWr-t \ n Rf,i 



R 



f,T-t+l 



1 + llj:_^j^^Yirj._^j^^jlrj._ 



(35) 



Corollary 3.2 shows that the solution of (23) at each time point T — t is similar to the 



solution of the common single-period optimization problem given in (33). The only difference 



is that the coefficient o^r-t is time varying. It depends on the future returns of the riskless 
asset. Leippold et al. (2004) suggested a similar solution to one presented in Corollary 3.2 by 
using the geometric approach in the case of independent asset returns. 

It appears that the expressions presented in Theorem 3.1 are quite complicate and cannot be 
evaluated for an arbitrary model of the asset returns. For this reason we suggest a reasonable 
approximation for the weights and evaluate its accuracy. 



First, the moments presented in (30) and (31) are approximated. Let ^ be a random variable 



with support [0, 1]. The application of the Sherman-Morrison formula to 1 — ST-t+2 leads to 

1 + (1 - Ol^T-t^l^T-t+lt^T-t+l 1 + (1 - O^T-W 



1 - 



1 + /j^T-t+l^T-t+lf^T-t+l 



1 + ST-t+2 



(36) 



where ST-t+2 = P^T-w^T-wt^T-W' Let r]T-t+2 = z whose support is [0,1] 

1 + ST-t+2 

as well. Then, it holds for the i-th component of that 



Er-t ( (1 — C^T-t+2)XT-t+l,, , . ^ 

= Et-I (^r]T-t+2^T-t+l,i^ = Et-I [r]T-t+2 {^T-t+l,i - t^T-t+l.i + t^T-t+l.i 
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Because rjT-t+2 < 1 it holds 



X 



< 



< ^/MSET-t+l^^ , 



where MSET-t+i,i = Et-i \ 0^T-t+i,i — f^T-t+i^iTj the conditional mean square prediction 
error calculated for flrj._^^^^.. If MSET-t+i,i is small, what should be expected if a good forecast 
of the process future values is performed, then the quantity ^ is well approximated by 



Similar results are obtained for (30). Here, it holds that 

Er-t (^(1 — ST-t-\-2)^T-t-\-l,i^T-t-\-l,j^ — Er-t (^VT-t-\-2^T-t-\-l,i^T-t-\-l,j^ 



(37) 



^iEr-t [r]T-t+20^T-t+i,j - Mt-^+ij)) + i^T-t+ijEr-t (vT-t+2{'^T-t+i,i - 



- ipT-t+l,ij - ^T-t+l,ij)ET-t (^T-t+2) + Er-t (vT-t+2{'^T-t+l,i - MT-^+l,i)(XT-t+l,j - Mt-^+Ij) 
+ flT_^^^^.ET-t {r]T-t+2{^T-t+l,j - MT-t+lj)) + f^T-t+ljEr-t (^^T-t+2(XT-t+l,i - /XT-t+l,i)) i 

where B^-t+i = ET-tO^r-t+i^T-t+i) = {bT-t+i,ij)ij=i,...k and S^-f+i = {crT-t+i,ij)ij=i,...k' 
Because < r]T-t+2 < 1 we get that 



En 



-t [pT-t+2{^T-t+l,i - MT-t+l,i)(XT-t+l,j - MT-t+lj)) - CFT-t+l,ijET-t (^T-t+2) 



WT-t+l,ij\ 

Er-t (\^T-t+i,i - MT-t+i,il) 



+ flT-t+l,iET-t {r]T-t+2{^T-t+l,j - At-^+Ij)) 

+ flT-t+l,jET-t {pT-t+2{^T-t+l,i - Mt-^+I,*)) 

+ \t^T-t+l,i\ET-t (l^T-t+lj - flT-t+l,j\) + IMT-^+IjI^'t 

< ^MSET-t+l,i^MSET-t+l,j + kT-t+l,ii I 

+ ||MT-^+l,il\/^^^T-t+l,, + \^lT-t+l^JWMSET-t+l^^ . 

Hence, if | and MSEx-t+i^i are relatively small for all j = 1, k we get the following 

approximation 



(38) 
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t^- 11 ^ ^ ^ A 1 + (1 - ET-t{r]T-t+2))sT-t+l r , ^ o 

l^maiiy, we note that tjt = I and rjT-t+i = z tor t > 2. 

Putting (37) and (38) together we obtain the fohowing approximation for the weights 




The approximation (39) can be used for arbitrary models of the asset returns for which 
\aT-t+i,ij\ and MSET-t+i,i are relatively small for all i^j = Moreover, this approx- 

imation becomes an exact one if the asset returns are independently distributed as it is shown 
in Corollary 3.2. 

Very interesting results can be obtained for the tangency portfolio as well (see, e.g. IngersoU 
(1987, p. 89), Britten- Jones (1999)). The weights of the tangency portfolio are derived at 
each time point T — t from ([34]) by assuming that the whole wealth is invested only into the 
risky assets, i.e. under the assumption w^_^l = 1. The results for the multi-period portfolio 
allocation problem with the tangency portfolio under the assumption of independence are given 
in Theorem 3.3. 

Theorem 3.3. Let = (X^^i,X^^2, • • • ,-^r,fc)^ r = 0, . . . be a sequence of the indepen- 
dently distributed vectors of k risky assets^ let rf^^ be the return of the riskless asset, and let 
E(X.r) = /x^ and ]/ar(X^) = where is positive definite. Then for all t = 1^ . . . the 
tangency portfolio weights for period T — t are given by 

^TP,T-t - 7^1 - -,,^-1 Z • l^Uj 

Theorem 3.3 shows that for the tangency portfolio the multi-period portfolio selection prob- 
lem is equivalent to the single-period allocation problem solved at each time point T—t. Because 
the tangency portfolio is, usually, considered as a market portfolio in the single-period alloca- 
tion problem (see, e.g., Britten- Jones (1999)), it can also be treated as a benchmark portfolio 
in the multi-period case, provided that the asset returns are independent. 



4 Empirical Illustration 

In this part of the paper we apply the results of Section 3 to real data. In the first example we 
deal with the multi-period portfolio selection under return predictability. In this subsection a 
model presented by Brandt and Santa-Clara (2006) is applied and approximative expressions to 
those presented in Theorem 3.1 are analyzed. In the second empirical illustration we consider 
an investor who attends to invest into the international portfolio consisting of five developed 
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stock markets, namely Belgium, Germany, Japan, the UK, and the USA. The approximative 
solutions of multi-period portfolio selection problem are obtained in this case as well. 



The investment strategy based on equation (39) will be called as the investment strategy 
based on the predictive loss approximation with a riskless asset and it is briefly denoted as the 
LAMPS strategy. It is noted that the LAMPS strategy coincides with the expressions given 
in Corollary 3.2 that are derived for the multi-period portfolio selection problem with riskless 
asset under the assumption that the asset returns are independent. Similar results are obtained 
in the case of the multi-period tangency portfolio (MTP). The application of the predictive loss 
approximation leads to the formulas presented in Theorem 3.2. 

As a flrst benchmark strategy of the multi-period portfolio selection we consider an investor 
who chooses the global minimum variance (GMV) portfolio which in the present situation is 
calculated by the time-invariant weights expressed as 

^GMV^t = -^,^-1-^ fol" t=l,...,r. (41) 

Because the same proportions of the wealth are invested into each stock at each time point we 
consider this portfolio as a myopic strategy. 

We also compare the results obtained by using the predictive loss approximation to the 
expression of the weights presented in Theorem 3.1 with two other benchmark strategies in the 
case with a riskless asset. The flrst strategy is known as a partial myopic strategy, i.e. the whole 
wealth is invested into the riskless asset. This strategy was suggested by Mossin (1968) who 
derived conditions under the utility function under which this strategy is optimal in the case 
of the portfolio consisting of one stock and the riskless asset. A further considered benchmark 
method is one suggested by Brandt and Santa Clara (2006) that is based on the assumption 
that 

w, = ezt , (42) 

where Zt is the vector of predictable variables taken at time point t. We refer to this portfolio 
strategy as the BSC strategy. In the empirical illustration, Brandt and Santa Clara (2006) 
argued that the solution of the multi-period portfolio selection problem based of this approach 
is very close to the one obtained from the simulation approach applied to the dynamic portfolio 
choice problem as suggested by Brandt et al. (2005) and van Binsbergen and Brandt (2007). 
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4.1 Multi-Period Portfolio Selection under Return Predictability 



In this section we deal with the multi-period portfolio selection problem assuming that the asset 
returns are predictable. This is one of the most commonly used approach applied for modeling 
the time series properties of the asset returns (see, e.g. Campbell and Viceira (2002), Brandt 
and Santa Clara (2006)). Here, we consider a special case of a model suggested by Brandt 
and Santa Clara (2006). They applied it to monthly returns from January 1945 to December 
2000 of one stock (r^) and one bond (r^). The term spread is used as a predictable variable Zt. 
Brandt and Santa Clara (2006, p. 2200) obtain the model 



(43) 
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(44) 



Note that the model (|43j) can be presented as a vector autoregressive process of order 1. Let 
Yj = (ln(l + rj%i), ln(l + r^^J, z^' then (lisj) becomes 



with Et ~ A/'(0, S). Thus, for Xj we get an autoregressive representation with 
Xt = Li> + L$Yt_i + Le^ = 1/ + $Yt_i + with L = [I 0] , 



(45) 



(46) 



where I is a 2 x 2 indentity matrix and is a 2 x 1 vector of zeros. Consequently, X^|J^t_i ^ 
A/'(/x^, S). We are interested in the conditional mean vector and in the conditional covariance 
matrix of given Tt-\. Note that 



iu, = ^(Xil Ji_i) = 1/ + *Yi_i, S = Var(Xi| j;_i) = LSL' , 



(47) 



Inserting (47) in (39) leads to the weights of the LAMPS strategy. In similar way the weights 
of the other multi-period portfolio strategies are calculated. 

We compare the performance of the above derived strategies with each other via an exten- 
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sive simulation study based on 10^ independent repetitions. Multi-period portfolio strategies 
are constructed for T G {6,12,18,24} and a G {0.833,0.909,0.937,0.952}. The values of a 
correspond to 7 G {5, 10, 15, 20} which are also used in Brandt and Santa Clara (2006). 



T ^--^ 


5 


10 


15 


20 


Method 


5 


0.5904 (0.0061) 
0.5847 (0.0065) 

0.5837 
0.5834 (0.0157) 


0.5474 (0.0016) 
0.5461 (0.0031) 

0.5456 
0.5459 (0.0039) 


0.5324 (0.0008) 
0.5317 (0.0017) 

0.5316 
0.5317 (0.0017) 


0.5245 (0.0004) 
0.5239 (0.0012) 

0.5241 
0.5242 (0.0010) 


LAMPS 
GMV 
Part. Myopic 
BSC 


12 


0.5949 (0.0046) 
0.5858 (0.0092) 

0.5839 
0.5822 (0.0172) 


0.5487 (0.0012) 
0.5460 (0.0035) 

0.5457 
0.5455 (0.0043) 


0.5330 (0.0006) 
0.5310 (0.0024) 

0.5316 
0.5316 (0.0019) 


0.5249 (0.0003) 
0.5229 (0.0021) 

0.5241 
0.5241 (0.0011) 


LAMPS 
GMV 
Part. Myopic 
BSC 


18 


0.5975 (0.0027) 
0.5867 (0.0106) 

0.5841 
0.5801 (0.0195) 


0.5493 (0.0007) 
0.5452 (0.0044) 

0.5458 
0.5450 (0.0049) 


0.5333 (0.0003) 
0.5298 (0.0035) 

0.5317 
0.5313 (0.0022) 


0.5250 (0.0002) 
0.5218 (0.0032) 

0.5242 
0.5239 (0.0012) 


LAMPS 
GMV 
Part. Myopic 
BSC 


24 


0.5989 (0.0014) 
0.5869 (0.0113) 

0.5843 
0.5757 (0.0240) 


0.5497 (0.0004) 
0.5441 (0.0054) 

0.5459 
0.5441 (0.0058) 


0.5335 (0.0002) 
0.5285 (0.0047) 

0.5318 
0.5309 (0.0027) 


0.5251 (8.5 ■ 10~^) 
0.5204 (0.0044) 

0.5242 
0.5237 (0.0015) 


LAMPS 
GMV 
Part. Myopic 
BSC 



Table 1: Medians and median absolute deviations (MAD, in parentheses) of the expected 
quadratic utility for the LAMPS, the multi-period GMV, the partial myopic, and the BSC 
strategies in the case of data from the example of Brandt and Santa Clara (2006). 



In Table 4.1 we present the median of the expected utility functions for the considered 
values of 7 and the investment horizon T. In each block the strategies order is the LAMPS, 
multi-period GMV, partial myopic, and BSC. In the parentheses the median absolute deviation 
(MAD) is given that is calculated as the median of the absolute deviation of the sample values 
from the median. The monthly data from January, 1945 to December, 2000 from the example 
of Brandt and Santa Clara (2006) are used that consist of one stock and one bond return at 
each time point the portfolio, while the term structure is used as a predictable variable for 
modeling time series properties of the return process. 

We observe a very good performance of the LAMPS strategy which is the best strategy 
for all T and 7. On the other hand it is not possible to provide a clear ranking between the 
GMV, partial myopic, and BSC strategies. While the GMV strategy is on the second place for 
smaller values of 7 and each T, the opposite results are observed for larger values of 7 where the 
partial myopic and the BSC strategies are ranked on the second and third places. Moreover, 
we obtained very small values of the MAD for the LAMPS strategy, while the largest values 
are attained for the GMV approach. 

In Figures 4.1-4.4 we study the problem in more details. Here, we present the empirical 
distribution functions of the expected quadratic utility functions obtained via the simulation 
study. If we compare several strategies based on the performance of their empirical distribution 
functions we should choose the strategy whose distribution function lies below the others. 
It follows from the fact that in this case for a fixed value of W the probability to obtain 
stochastically larger value is the largest one, i.e. with the higher probability larger values of the 
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expected utility are attained. Because the partial myopic strategy is based on investing into 
a riskless asset it is a deterministic strategy and it's distribution function is a step function. 
The position of the given index is used for comparing the dynamic portfolio strategies. This 
criterion is based on the maximum probability of exceeding the utility of the partial myopic 
portfolio. 

In Figure 4.1 we observe a very good performance of the LAMPS and BSC strategies, 
especially for smaller values of 7 and T. If 7 increases then the probability of overperforming 
the partial myopic strategy decreases. It is always small for the GMV approach. For this 
method it is less than 40% for larger values of 7 and T, while it is very large for the LAMPS 
strategy (about 80%). 

4.2 Multi-Period International Optimal Portfolio 

We consider an investor who invests into an international portfolio consisting of the capital 
markets indices of five developed stock market, namely Belgium, Germany, Japan, the UK, 
and the USA. The weekly data of the MSCI (Morgan Stanley Capital International) indices 
for for the equity markets returns are considered from the 4th of January, 2002 to the 4th of 
December, 2009 and the return series for each index are calculated. To the return series we fit 
the VAR(l) process defined by 



Xt = + ^Yt-i + St with St ^ uM{0, S^) . 
We obtain 



(48) 



u — 
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-04 




0.2011 
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0.455 


1.20e 


-03 




0.3139 


-0.1231 


-0 


00191 


-0 


511 


0.434 


6.74e 


-04 




0.0487 


0.0888 


-0 


12131 


-0 


224 


0.343 


5.54e 


-04 




0.1829 


-0.0889 





00988 


-0 


441 


0.382 


2.79e 


-05 




0.0766 


-0.0643 


-0 


03049 


-0 


114 


0.133 



and (49) 



0.0013085186 0.0010544496 0.0004365753 0.0009120373 0.0006781289 

0.0010544496 0.0013833540 0.0005648237 0.0010218539 0.0008332314 

0.0004365753 0.0005648237 0.0007994341 0.0004733366 0.0003667012 

0.0009120373 0.0010218539 0.0004733366 0.0010176793 0.0006927251 

0.0006781289 0.0008332314 0.0003667012 0.0006927251 0.0007242233 

It is noted that the last column of the matrix $ has the largest values. It shows that the 
influence of the US market on the return indices is larger than those of the domestic ones. 
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Next, calculate the weights of the three multi-period portfolio strategies (LAMPS, GMV, 
and partial myopic) as already described in Section 4.1. We exclude the approach of Brandt 



and Santa Clara in this study since there is no predictable variable within the model (49). We 
choose 7 G {5, 10, 15, 20} and T G {4, 8, 12, 16}. 

In Table 4.2 we present the medians and the mean absolute deviations for the expected 
utilities obtained via a simulation study with 10^ independent repetitions of the autoregressive 



process (48) with parameters (49). The results are similar to those of Section 4.1. We observe 
a very good performance of the LAMPS strategy which turns out to be the best one. On the 
second place the partial myopic strategy is ranked followed by the multi-period GMV portfolio. 



^\ 7 
T ^\ 


5 


10 


15 


20 


Method 


4 


0.5988 (0.0014) 
0.5832 (0.0052) 
0.5836 


0.5497 (0.0004) 
0.5454 (0.0028) 
0.5456 


0.5334 (0.0002) 
0.5314 (0.0017) 
0.5316 


0.5251 (9.5 • 10-^) 
0.5239 (0.0012) 
0.5240 


LAMPS 
GMV 
Part. Myopic 


8 


0.6002 (5.9 • 10-^) 
0.5829 (0.0074) 
0.5838 


0.5500 (1.5 • 10-^) 
0.5452 (0.0037) 
0.5456 


0.5336 (6.6 • 10-*^) 
0.5311 (0.0022) 
0.5316 


0.5252 (3.4- 10-*^) 
0.5234 (0.0016) 
0.5241 


LAMPS 
GMV 
Part. Myopic 


12 


0.6002 (2.1 • 10-^) 
0.5826 (0.0090) 
0.5839 


0.5501 (5.2 • lO-'O 
0.5449 (0.0043) 
0.5457 


0.5336 (2.4- lO-'O 
0.5306 (0.0027) 
0.5316 


0.5252 (1.5 • lO-'O 
0.5229 (0.0021) 
0.5241 


LAMPS 
GMV 
Part. Myopic 
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0.6002 (1.4- lO-'O 
0.5824 (0.0102) 
0.5840 


0.5501 (8.4- 10-^) 
0.5445 (0.0048) 
0.5458 


0.5336 (6.7- 10-^) 
0.5302 (0.0031) 
0.5317 


0.5252 (5.7- 10-^) 
0.5224 (0.0026) 
0.5241 


LAMPS 
GMV 
Part. Myopic 



Table 2: Medians and median absolute deviations (MAD, in parentheses) of the expected 
quadratic utility for the LAMPS, the multi-period GMV, and the partial myopic strategies in 
case of weekly MSCI data for the period from the 4th of January, 2002 to the 4th of December, 
2009. 



More pronounced results are presented in Figure 4.5. Here, we observe a very good per- 
formance for the LAMPS strategy for all 7 in the case of T = 12 weeks. Both strategies 
overperform the partial myopic strategy with probability of almost equal to one if 7 = 5. For 
larger values of 7 a similar behavior is present for the LAMPS approach. The multi-period 
GMV portfolio performs much worse. For 7 = 5 and 7 = 10 the probability of getting higher 
values of the expected utility is about 50%, while it is less than 40% for 7 = 5 and 7 = 10. 
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Next we want to present a further example. Now the return of the US market index is used 
as a predictable variable and an investment into four other capital market indexes is considered. 
This example is motivated by economic theory - the influence of the US market may be larger 
than of the domestic one, and it is justifled in the structure of the matrix $ (cf. ([49|)). In this 
example all of the above described strategies are analyzed. 

In Table 4.3 we present the medians and the MADs of the expected utilities. A very good 
performance of the LAMPS multi-period portfolio strategy is observed. This approach is ranked 
on the flrst place. A much worse performance can be observed for the partial myopic strategy, 
for the multi-period GMV portfolio, and for the BSC method which can be ranked on the 
second, third, and fourth places. 



^\ 7 
T ^\ 


5 


10 


15 


20 


Method 


4 


0.5981 (0.0021) 
0.5835 (0.0059) 

0.5836 
0.4621 (0.1270) 


0.5495 (0.0006) 
0.5455 (0.0031) 

0.5456 
0.5172 (0.0306) 


0.5333 (0.0003) 
0.5315 (0.0018) 

0.5316 
0.5187 (0.0139) 


0.5251 (0.0001) 
0.5238 (0.0013) 

0.5240 
0.5170 (0.0077) 


LAMPS 
GMV 
Part. Myopic 
BSC 


8 


0.6001 (0.0002) 
0.5835 (0.0083) 

0.5838 
0.5309 (0.0634) 


0.5500 (3.97- 10-^) 
0.5454 (0.0039) 

0.5456 
0.5338 (0.0151) 


0.5336 (1.9 • 10-^) 
0.5310 (0.0023) 

0.5316 
0.5264 (0.0067) 


0.5252 (10-^) 
0.5233 (0.0018) 

0.5241 
0.5212 (0.0037) 


LAMPS 
GMV 
Part. Myopic 
BSC 


12 


0.6002 (9.8 • 10-^) 
0.5835 (0.0100) 

0.5839 
0.5370 (0.0578) 


0.5501 (2.4- 10-^) 
0.5451 (0.0044) 

0.5457 
0.5349 (0.0140) 


0.5336 (1.1 • 10-^) 
0.5304 (0.0029) 

0.5316 
0.5271 (0.0060) 


0.5252 (5.8 • lO-'O 
0.5226 (0.0024) 

0.5241 
0.5215 (0.0034) 


LAMPS 
GMV 
Part. Myopic 
BSC 


16 


0.6002 (5.7 10-0 
0.5836 (0.0113) 

0.5840 
0.5355 (0.0593) 


0.5501 (1.4- lO-'O 
0.5445 (0.0049) 

0.5458 
0.5344 (0.0144) 


0.5336 (6.3 • 10-^) 
0.5298 (0.0035) 

0.5317 
0.5269 (0.0063) 


0.5252 (4- 10-^) 
0.5219 (0.0030) 

0.5241 
0.5214 (0.0036) 


LAMPS 
GMV 
Part. Myopic 
BSC 



Table 3: Medians and median absolute deviations (MAD, in parentheses) of the expected 
quadratic utility for the LAMPS, the multi-period GMV, the partial myopic, and the BSC 
strategies in the case of weekly MSCI data for the period from the 4th of January, 2002 to the 
4th of December, 2009. The investment into the four countries (Belgium, Germany, Japan, and 
the UK) are considered, while the returns of the USA market are used as a predictable variable. 



The same ranking is also observed in Figure 4.6 where T = 12 is used. It has to be 
emphasized that the LAMPS strategy overperforms the partial myopic strategy with probability 
almost 1 for all of the considered values of 7, while the multi-period GMV portfolio and the 
BSC method are really bad. 
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5 Summary 



Although in 1959 Markowitz has already formulated the multi-period asset allocation problem, 
no closed-form solution is available in literature up to now. Brandt and Santa Clara (2006) 
provided an analytical solution by imposing some constraints on the structure of the portfolio 
weights and, thus, by transforming the multi-period portfolio selection problem with time- 
dependent weights into a more simpler one where the matrix of time-independent constants 
has to be calculated at the beginning of the investment horizon. Moreover, note that the 
approach may lead to a local maximum and not a global one. 

In the present paper, we derive a closed- form solution of the multi-period portfolio selection 
problem with and without a riskless asset which is obtained under weak conditions on the 
process of the asset returns. The expressions of the optimal portfolio weights are based only 
on the conditional mean vectors and the conditional covariance matrices. Moreover, under 
the additional assumption of independence they are similar to the solutions of the single-period 
optimization problems that are performed at each time point. We prove that only the coefficient 
of the shape parameter depends on the dynamics of the asset returns in the case without a 
riskless asset. In the presence of a riskless asset the solutions of the multi-period problem 
and of the sequence of simple-period problems are proportional where only the coefficients of 
proportionality depend on the process dynamics. A very remarkable result is obtained for the 
portfolio selection problem based on the tangency portfolio. Assuming independent returns it 
is proved that the solution of the multi-period problem and the solutions of the simple-period 
problems are the same. 

The derived multi-period portfolio choice strategies are compared with the existing methods 
proposed in literature for two real data sets. The first data are taken from the paper of Brandt 
and Santa Clara (2006), while the returns of five developed stock market indices are used in the 
second example. We observe a very good performance of the predictive loss approximation of 
the multi-period strategy in case with a riskless asset which is always ranked on the first place. 

The obtained results can be further extended by allowing the intermediate consumptions. 
This can be done without any large effort by slightly modifying the weights of Theorems 2.1, 
3.1 and 3.2. 

6 Appendix 

In this section the proofs of the theorems are given. 
Proof of Theorem 2.1: 
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First, we note that the expression of the optimal weights at period T — 1 is given in (11). The 
rest of the proof is done by using the mathematical induction on the expressions of the portfolio 
weights and the value function. Let 



A, = E,_i[V,+i±,±[] for i = 1, . . . , T - 1 and At = 'Et + P^tP^t • 



(50) 



Moreover, let 



Mi 



Mt 



for i = T 



Ei_^[Ri+iXi] for ^ = l,...,^-l, 



(51) 



Ri = ^ , Vi = .-^-1 and Si = fx* 'Qip,* with 

1 -/x^- 1 1 -/I, 1 



= A"' - ^r^3T^ for t = i,...,T-i 



Note, that 



1U = 0^ Q^1 = 0, Q.A,Q, = Q,. 



Let be the optimal portfolio weights calculated for period T — 1 given in (11). First, 

we calculate the value function for period T — 2. It holds that 



V{T-2,Wt-2.Tt-2 

- W7-_i:w^_-|^l=l 



max Et-2 

W7-_2:w^_2l = l 



max [Wt-i^t-iPt ~ —W^_iW^rj._^ATV^T-i 



max Et-2 

W7-_2:w^_2l = l 



max Et-2 

W7-_2:w^_2l = l 



Wt-i^t'-iPt ~ -:^W^_iW^'_iAt^t-i 

, ^T^l 1 - \ 

Wt-1 I + -777 QtAt At 



VA-^l aWr-i 



-Wr 



T-l 



A-^l 



A-^l 



oWt- 



QtPt 



( 



max Et-2 

W7-_2:w^_2l = l 



Wt- 



\ 



\ 



I'A'H ~ 2 I'A'H 2 I'A'H 
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Using the definitions of i?^, 1^ and Si we obtain 



y(T-2,WT-2,-FT-2) 



max Et-2 

W7-_2:w^_2l = l 

max Et-2 

W7-_2:w^_2l = l 



1 



Wt-iRt + - ^W^_^Vri 



2a 



T 



2a 2 



max 



a 

WT-2^T-2i^T-l + F{st) - -^Wt_2 (w^_2 A^-l WT-2) 



where 



F{St) = —Et-2[st] 
za 



F{st) does not depend on wt-2- 

The last expression is similar to the value function at period T — 1 (cf. Hence, 
optimal weights w^_2 are given by 



w 



T-2 



H — — — Qr-iAr-i 'with Qt-i = ^r-i ~ 



l'A^\l aWT-2 



As a result, the following expressions are the basis of the induction 



V{T-2,Wt-2,J't-2) 



w 



T-2 



max 

W7-_2:w^_2l = l - 

a 



WT-2^T-2i^T-l + F{St) 



Wt-2^'t-2^T-1^T-2 



A-\l 1 

H 77T Qt-IMt-I 



VA-'_{\. aWT-2 



with F{st) as defined in (52). 

In the induction hypothesis we assume that the statement holds for t = i.e., 



V{T-n,WT-n,J'T-n) 



max 

WT-n:W^_^l = l 

a^ 



WT-n'^T-nf^T-n+l 



W 



T-n 



1 1 



Qr -n+lMT-n+l ' 



where 



1 / 

F(st, St-1, Sr-n+2) = ^ I Et-2[st] + ^m-2[ 



m=T-n+2 
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Note that the last quantity does not depend on w^-n- 

In the inductive step we prove that the last identities also hold for t — n + 1. It is sufficient 
to derive the value function for period T — (n + 1) which is given by 



V{T - {n + l),WT-^r.+l),J'T-(n+l)) 



^T-(n+l) {WT-n^T-nP'T-n+1 + ^(^T, ^T-l, • • • , 5T-n+2) 



max ET-fn+i)\WT-n 
^T-(n+l)l = l ^ 



VA, 



T-n+l-" 



2 VA 



T-n+l-" 



2 l^A 



T-n+1' 



\ —F^T-n + l 



= ST-n+l 



Hence, we obtain 



V{T - (n + 1), Wr-(n+l), ^T-(n+l)) 



max 



^T-(n+l) 



1 = 1 



(n+1) 



max 

^T-(n+l) 



2^T-(n+l)W; 



T-(n+l)^T-nWT-(n+l)^ 



where 



It is the desired form of the value function at period T — (n + 1). Because this expression is 
similar to the value function at period T — we get the following formula for the weights at 
period T — (n + 1) 



w 



T-(n+l) 



1 



VA-\1 aW, 



1 



T-(n+l) 



-Qr-nP^T-n 



The theorem is proved. 



For proving Corollary 2.2 we use the result of Proposition 6.1. 



Proposition 6.1. Let X &e a random vector with mean /x and positive definite covariance 
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matrix S. Let A = S + /i/i' with /i = /x + 1. // 

w ^ — h fi"-^Q/i TOi/i Q A"-^ ^ (54) 



w ^ + a"^Q/i TOt/i Q S"^ ^ (55) 



and 



~ (i + A's-V)i'S-4-(i'S-iA)' " T+l ' ^ ^ 

where Rqmv ^ y^-^ ^ ^ fi'^lfi li'CllJi. 
Proof of Proposition 6.1: 



From (55) we obtain 



where 



Ci = ^-C2^^^^ and C2 = a-\ (58) 



In order to prove the proposition we need to show that (54) can be expressed in the same 



way. The apphcation of the Sherman-Morrison formula (HarviUe (1997, Theorem 18.2.8)), i.e.. 



A-' = (S+A/i')- = S-'-^-^" 

1 + /X 2] /X 



leads to 



w = {l-Ka-^)- /^^T'^. (59) 

^ ^(l + /i'S-V)l'S-4-(l'S-^/i)2 ^ ' 

/ /i'S"^l(l - Kd-l) ^y;-l~ 

^ V (1 + - (I'S-^A)' ^ 1 + / ^' 

where 

1 + /Li E IX 
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From the structure of (57) and (59) we get 



a 



Co 



(1 + /i's-V)i'S-4 - (l'S-^/i)2 



+ 



a-i - 1 - R 



GMV 



For proving the proposition we only need to show the equahty of the coefficients in front of 



in (|57|) and (|59|). It holds that 



Ci 



a 



Co 



1 + A'E-i/i (1 + /i'E-V)!'^:-^ - (i'5]-i/i)2 J I's-^i 



(l-Ka 
(1 - Ka-^ 



+ 



i'S-4 ((1 + - (i'S-i/i)2) 

(1 + 

(1 + /i's-i/i)i'S-4 - (i'S-i/i)2 



The last identity completes the proof. 



Proof of Corollary 2.2: 

Under the assumption of independence it holds that 

Et + AtAt for ^ = 1 .^-.x 

At-«+i = <( . (61) 

yT-t+2{'^T-t+i + Ar-i+iAr-t+i) fo'^ t = 2, . . . ,T , 
and 

{jUt, for t = 1 

(62) 
RT-t+2iiT-t+i for t = 2,...,T 

Let i?T+i = ^+1 = 1- Then, 

(VT-i+2(^T-m + AT-t+lA0)~^l , 1 ^ ~ ^ 

WT-t = ^TTTT? 7^^ — r ^, Y^ + ~iri ^T-t+l{nT-t+2fiT-t+l) 

i- [yT-t+2[^T-t+i + ^iT-t+l^^T-t+l)) i aWT-t 

^ (^r-t+1 + jj'T-t+lP'T-t+l)'^'^ , RT-t+2 A 
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with 



and 



where the last identity foUows from the definition of i?T-t+2 and VT-t+2 given in (16). 

The rest of the proof foUows from Proposition 6.1 if S is replaced by S^-t+i, A by jjiT-t+i 
and 



a 



aWT-tVT-t+2 aWr-t \ . ^ ^ aWr-. v • .^o 

\i=l —t-\-2 / \i=l —t-\-2 



T 



where 



a, = l (E, + /x,/xj /X, = . 

(1 + RgMV^iY + (1 + SijVcMV.i 

The last expression is obtained by applying the Sherman-Morrison formula. At last, we recall 



Qr-t+il = and get (18). Thus the corollary is proved. 



Proof of Theorem 3.1: 



The expression of the optimal weights at period T — 1 is given in (28). The rest of the theorem's 
statement is proved by using the mathematical induction on the expressions of the portfolio 
weights and the value function. We use similar notations as in the proof of Theorem 1. Let 
Ai = Ei_i[{l - 5i+i)X^X.] for i = 1, . . . , T - 1 and At = St - fJ^rf^'r^ 

J I^T for i = T 

l^i = \ - (63) 

^z-i[(l-5,+i)X,] for z = l,...,r-l, 



and Si = /i* ^A- /i* for i = 2, T. 

Let Wt^i be the optimal portfolio weight calculated at period T — 1 in the case of a riskless 
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asset as expressed in (28). First, we calculate the value function at period T — 2. It holds that 



V{T -2,Wt-2.J^t-2) 

1 



max 



aax Et-2 \Wt-1 ( Rf,T + ( —r- %,T ) t^'T^T^ I^T 

T-2 L \ \aWT—i / 



2 ^"^ WaWr-i 



\aWT- 



—1 ^ 2 



2i? 



f,T I — %,T ) flT^T^t^T 



y aWT—1 

Using the definition of st we obtain 



V{T -2,Wt-2^J^t-2) = mdixET-2 Wr-iRf t{1 - st) 

wt-2 L ' a 



St 



2 I Icel^T-i 



aWr-i 



Rf T ] St 



max £^T-2 

Wt" — 2 



max 

Wt" — 2 



2a 2 ^ 



^T-2^/,T (w^_2At-1Wt-2 + £^T-2[1 " ^t]^/,T-1 + '^Rf,T-l'^T-2f^T-l 



Et-2[st] ' 

2a 



The last expression is similar to the value function at the period T— 1. Hence, it is maximized 
on the weights w^_2 expressed as 



w 



T-2 



T-2 



(64) 



Hence, the basis of induction are the following expressions 



V{T - 2, Wt-2, = max VI/t-2^/,t (&t^/,t-i + w^_2/i^_i) 



2 
w 



Wt_iR)^T (w^_2At-iWt-2 + bTRf^T-1 + 2i?/,T-i w^_2/i^_i) + F(5t) 

1 



T-2 



T-2 



with F{st) = ^^f^ and 6t = ^t-2[1 - ^r]. 
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In the induction hypothesis we assume that the statement holds for t — i.e., 

V{T -n,WT-n,^T-n) = , max Wr-n TT Rf,i {bT-n+2Rf,T-n+l + "Vv'j.-nfJ'T-n+l) 

— > ^ \i=T-n+2 J 

\i=T-n+2 J 
+2) , 

^T-n = ^^^^^^^ ^ n -%,T-n+l^ A^i^^l/i,5._^^l . 

/ T-l T-1 \ 

with F{st, . . . , 5T-n+2) = ^ ^t-2[^t] + E 11 &2^m-2[5m] and bi = £;i-2[l - 5^]. 

V m=T-n+2 i=m J 

In the inductive step we prove that the last identities also hold for t = n + 1. It is sufficient 
to derive the value function for period T — (n + 1) which is given by 

V(T - (n + 1), + -^T-Cn+l)) 

= ^ max Ej._^^_^^^^WT-n [ Yl ^f,i ] (^T-n + 2Rf,T-n + l + ^T-nP-T-n + l) 
T-(n+l) \i = T-n + 2 J 

- ^W^_^ I Yl (wT-n^T-n+lw5._^ + 6T-n + 2J?/,T-n + l +2J?/,T-n + lW?.i^/iT-n + l) 

^ \i = T-n + 2 I 

F(st, . . . , ST-n + 2)] 



AT-n + l^T-n + l'^T-n+l) 

T-n \i = T-n + 2 J J 



T 

Using the definition of Si and denoting ^ — Yl we receive 

V(T - (n + 1), TyT-(n+l),^T-(n+l)) 
= max £;T-(n+l) h^T-n%,T-n+l^^T-n+2(l " ST-n+l) + ^ "^^^ ^T-n+1 

+ F(ST, . . . ,ST-n+2) - -VK|._^6T-n+2^^ %,T-n+l^ + ^/,T-n+l 

T-(n+l) 

max m^T-n^%,T-n + l^T-n+2(l - ST-n + l) 



+ 2i?/^T-n + l I - Rf,T-n + l ) ^T-n+l 



{wT-(n+l)} 



= max WT-(n+l)^%,T-n+l^T-n+2 (^T-(n+l)(l - ^T-n+l)%,T-n + W^_/ ^-^n/x^ 



-Vt^T_(n + l)(C%,T-n+l)^^T-n+2 (w^_ (^^-^^^ At-tz W^- (n+1) +^T-(n+l)(l " ^T-n + l )^/,T-n 
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where F{st, • • • , sVn+i) = F{st, • • • , sVn+2) + "-"^^^ &T-n+2. 

It is a desired form of the value function at period T — (n + 1). Because this expression is 
similar to the value function at period T — n, we get the following formula for the weights at 
period T — (n + 1) 

- f (^-^/,T-n+l)~^ ^ \ ^-1 . * 



n 1 



T 

Substituting ^ — Y\ Rf,i leads to the expression given in the statement of Theorem 3.1. 

i=T-n-\-2 

The theorem is proved. 



Proposition 6.2. Let X 6e a random vector with mean /jL and positive definite covariance 
matrix S. Let A — J2 + fifi and fi — ii — r/1. // 

w = 7"-^A pi (65) 

then 

w = 7"'S-^/x m^/i = , ^ 1. . (66) 

1 + /X S /X 

Proof of Proposition 6.2: 

The application of the Sherman-Morrison formula, i.e., 

A-' = (s + MM)- = s-'- f3^y' 

1 + /U S 

leads to 

what completes the proof of the proposition. 



Proof of Corollary 3.2: 

Under the assumption of independence 

Sr + AtAt for t^l 

Ar-t+i = { (67) 
(1 - ST-t+2)(Sr-t+i + (iT-t+ilJ'T-t+i) for t^2,...,T 



36 



and 




At 

(1 — ST-t+2)P'T-t+l 



for t^l 



for t = 2,...,T 



(68) 



Then the statement of the corollary follows from Proposition 2 if S is replaced by S^-t+i 
and /i by /Xj^.^,^!, and 



Proof of Theorem 3.3: 

The results of Theorem 3.3 follow Theorem 3.1 and the application of the Sherman-Morrison 
formula. 
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